Atom-dimer scattering for confined ultracold fermion gases 
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We solve the three-body problem of a quasi-one-dimensional ultracold Fermi gas with parabolic 
confinement length a± and 3D scattering length a. On the two-body level, there is a Feshbach-type 
resonance at a±/a m 1.46, and a dimer state for arbitrary a±/a. The three-body problem is shown to 
be universal, and described by the atom-dimer scattering length a a d and a range parameter b a d- In 
the dimer limit a±/a S> 1, we find a repulsive zero-range atom-dimer interaction. For a±/a <C — 1, 
however, the potential has long range, with a a d > and b a d 3> ow- There is no trimer state, and 
despite a a d — at a± /a ~ 2.6, there is no resonance enhancement of the interaction. 

PACS numbers: 03.75.Ss, 05.30.Fk, 03.65.Nk 



The recent experimental observation of the formation 
of dimers (molecules) in ultracold binary Fermi gases 
|l| has sparked intense excitement and activity among 
both atomic and condensed-matter physicists. Exper- 
iments are now able to probe in detail the BEC-BCS 
crossover regime by using magnetic-field tuned Feshbach 
resonances [3- E3- DJ ■ where the 3D scattering length a 
describing the s-wave interaction strength among differ- 
ent fermion species can be tuned almost at will. This 
crossover remains a difficult and long-standing challenge 
to theory, as there is no small parameter in the prob- 
lem |5j. Here we discuss the related but simpler prob- 
lem of an ultracold two-species (T,i) Fermi gas con- 
fined in a ID trap potential. On the two-body level, 
there is always a bound state (dimer), and one has a 
confinement-induced resonance (CIR) in the ID atom- 
atom scattering length a aa [|| Q , similar to a Feshbach 
resonance. Such a scenario appears to be experimen- 
tally feasible 0, 01 H S IE EI and could reveal inter- 
esting new physics. In this paper, we analytically solve 
the three-body problem for confined ultracold fermions, 
and compute the atom-dimer scattering length a a d and 
the potential range b ac i [see Eq. (|11H below]. Our results 
also determine parameters entering models that may be 
solved exactly by powerful many-body techniques in ID 
[13, see also Refs. [HII2- 

For simplicity, we assume the two fermion species to 
have the same mass mo- Under the harmonic transverse 
confinement potential U c {r) = ^moLu^(x 2 + y 2 ), with 
associated lengthscale a±_ = (2/i/mow^) 1 / 2 , there is a 
two-body bound state (dimer) with dimensionlcss bind- 
ing energy 



n B = (huix - E b )/2Huj ± 
determined by the condition 



C(l/2, fi B ) + ax/a = 0. 



(1) 



(2) 



Since the zeta function £(1/2, fi) is monotonic in Q, 
there is exactly one bound state for any given a±/a, 
although the 3D problem has a bound state only for 



a > 0. For a±/a — > — oo, the 'BCS limit' is reached, 
where Qb — (a/a±) 2 <C 1, and the dimer size is large, 
of the order aj_/|a|. In the tightly bound 'dimer limit', 
aj_/a — ► +oo, the dimer size is small, of the order a, and 
£Ib — (a±/2a) 2 3> 1. The analogue of the Feshbach reso- 
nance is then realized by the CIR. Solving the two-body 
scattering problem with just one open channel, the ID 
scattering length is 



2 



ax 
a 



C 



C 



-C(l/2) ~ 1.4603. (3) 



At low energies, this implies that one can use the ID 
atom-atom interaction potential V aa (z,z') = g aa S(z — 
z') with g aa = — 2H 2 /moa aa . The CIR then occurs for 
&aa = 0, corresponding to Qb = 1, and can be reached 
by tuning a± or a. 

A natural question then concerns the scattering prop- 
erties of the atom-dimer system, for instance, the scatter- 
ing length a a( i- This problem has recently been solved for 
the unconfmed 3D case by Petrov [ijj], where a a d ~ 1.2a. 
In the confined geometry, it is then interesting to ask (i) 
whether the scattering length a a d also shows CIR-relatcd 
resonant behavior, (ii) whether a trimer state could be 
possible in the confined geometry, and (iii) whether a uni- 
versal description in terms of simple two-body physics is 
always applicable. For the unconfined case, Efimov [l4j 
showed that this three-body problem is universal. More- 
over, there is no bound trimer state. Below we shall 
answer these questions in detail. 

We study the three-body problem (TTi) with two iden- 
tical fermions, and denote by xi (x2,3) the position of the 
I (the two |) particles. Next we perform an orthogonal 
transformation to variables (x, y, z) in order to decouple 
the center-of-mass coordinate z [jjj. Since the confine- 
ment is harmonic, U c remains diagonal in the positions. 
The three-body problem then reduces to 



m 



V 2 , + U C (X) - E) *(X) = -J2 V(t±)*(X), 



where X = (x, y) is a six-dimensional vector, x 



(4) 



(2xi — x 2 — x 3 )/2 and y = x 3 — x 2 . With these def- 
initions, the distances between the | particle and each 
1 particle are r± = \/3x/2 ± y/2. We adopt the pseu- 
dopotential approximation for the 3D interaction, V(r) = 
(4:Trh 2 a/mo)5(r)J^;(r-), which allows to incorporate inter- 
actions via boundary conditions imposed for vanishing 
distances between f and 1 atoms. For r± — > 0, this im- 
plies the singular behavior 



*(X) 



./(rj 



4irr± 



-(l-r±/o), 



(5) 



where the rx,± = x/2=F\/3y/2 are orthogonal to r±. We 
consider E < and write 

E = -2n B hui ± + h 2 k 2 /m , 

where the relative momentum hk of the atom-dimer com- 
plex is sent to zero later. 

The boundary conditions J5J) allow to write Eq. J3J in 
the form 

± J 

(6) 
where the two-particle Green function is 



G^ ) (ri,r 2 ,r 3 ,r 4 ) = ^ 



E Xl + -Ea 2 - E 



Here ^a denotes the eigenfunctions to the single-particle 
problem for cigenencrgy E\ . The quantum numbers A in- 
clude the ID momentum k, the (integer) angular momen- 
tum m, and the radial quantum number n = 0, 1, 2 . . ., 
which gives E x = hui±(2n + \m\ + 1) + h 2 k 2 /mn and 
0a = e 171 ^ R nm {p)e ikz with radial functions R nm [if. Us- 

(2) 

ing the fact that G E and the integration measure are 
both invariant under orthogonal transformations, we find 
from Eq. (jfJJ) 



m *(r,r i )=y ( ir' ± /(rl)[Gg ) ( 



r,r i ;0,r , 



-Gf{r/2 + V3r ± /2, v/3r/2 - r x /2; 0, r^)], 

where r = r_ and rx = rx,-- Next we implement the 
r — ► limit according to Eq. JSJ to obtain a closed equa- 
tion for /(rx). This limit can be directly taken for the 
non-singular second term in the integral above, while the 
first term contains the singular behavior necessary from 
Eq. iJSJ. Once this singular behavior is removed, one ob- 
tains a regular integral equation for /(rx) |13l il5l Il6j. 
It is convenient to transform from real space into the 
complete basis {ipx}, implying 



X' 

where we use [see also Eq. (J2J] 

C(Q) ---- ((l/2,0)-C(l/2,fi fl ), 

(a ± k/2) 2 +E x /2fmjj 



(7) 



and the matrix A x>X i is given by 

4ira± 



m 



dr ± iP* x (r ± )*l> x ,(-r ± /2)G E - Ex {V3T ± /2,0). 



Using the integral representation of Ref. [llj for the two- 
body Green function Ce(f, r'), A\,y can be evaluated 
explicitly Before analyzing Eq. J2J further, however, it 
is useful to perform a rescaling. So far, / has been taken 
as a function of rx = (p, z). However, for the asymptotic 
solution consisting of a dimer and a free atom, the atom- 
dimer distance d does not coincide with rx- With the 
dimer wavefunction < &o( r ) [6j, we expect 

*(r,rx) a $ (r)x(d), d = ( Xl +x 3 )/2-x 2 = r + r . 

In the r — > limit, "f ~ x(v / 3rx/2)/47rr, which estab- 
lishes a connection between d and rx from Eq. J5J) . After 
the rescaling rx — ► V3rx/2, / coincides with the scatter- 
ing solution x- Therefore, from now on, all wavevectors 
are rescaled by the factor 2/\/3- 

Let us then proceed by projecting the integral equation 
to the lowest transverse state (n = m = 0). Explicit 
calculations [lTj show that the higher states are negligible 
in the BCS limit and affect a a d only slightly in the dimer 
limit, see below. Taking into account the above rescaling, 
noting that only m — modes have nonzero overlap with 
the lowest state we find 



+ OG 



C{Q k )f k = 

n k = n 

with the kernel given by 



.4 



k,k' 



E 



dk ' A f 

~Z — -rt-k.k' Jk' ! 
Z7T 

3aj(fc 2 - k 2 ) 
16 



4-p 



(8) 



^ P + tt B + (ax/2) 2 h3fc 2 /4 + k 2 + k' 2 + kk'} ' 

i-1 (9) 

Following Ref. [llj, we now make an ansatz for the 

solution of the integral equation, 



f(k) = 2nS(k-k)+if(k,k)Y^ 



1 



k ± k + i0+ 



(10) 



with a regular function f(k,k). This ansatz gives the 
expected scattering state after Fourier transforming to 
real space, 



' k I - 



ft) 



Q, 



f(z) = e lkz + / (sgn(z)k, k) e 

In the low-energy limit k, k — > 0, on general grounds |18j . 
the expansion 

f(k, k)=-l + ikb ad + ika ad + 0{k 2 , P, kk) (11) 

applies, where a ad is the atom-dimer scattering length. 
From the analysis of model potentials, the length b ac [ is 



linked to the range of the effective ID atom-dimer po- 
tential. In particular, b a d — » if a zero-range (5-potential 
can be used for the effective ID atom-dimer scattering at 
low energy scales. 

Inserting the ansatz (|10|l into Eq. (JSJ, we obtain 



C(il k ) 



k 2 -k 2 



2ikf{k, k) - iV 



± 



+°° dk' A 



k-.k-' 



2tt k±k 



-Jik',1) 



f(k,k)A k - k + f(-k,k)A k - k 



A h 



k.ky 



(12) 



where V denotes a principal value integration. Finally, 
the analysis can be simplified considerably by letting 
k — > 0, i.e., by expanding Eq. I|12[) in k and keeping only 
the lowest order. At that stage, we switch to dimension- 
less momenta u,v! by writing k = (2\ZTTg/aj_)u. Some 
algebra gives, with the weakly O^-dependent functions 



G{u,u') 
H(u) 



4-p 



p=0 



■p/Qi 



A~p U 



p=0 y 



■p/SIb) 



-2 ' 



the following integral equation for h(u) = f(u, 0): 

n+ °° du' 

[G(u,u')h(u') - G(u,0)h(0)} 

3u 2 ' 




4 



h(u) 



(13) 



G(u,0) + iH(u). 



Note that the real (imaginary) part of h(u) is even (odd) 
in u. The scattering length a ac i finally follows from the 
real part of Eq. (|13|1 and the condition h(0) = — 1, see 
Eq. (jllfl . while b ac i can be extracted from the imaginary 
part of Eq. I|13|) . The integral equation i|13|) shows in 
particular that a a d/a±_ depends only on Ob, and hence 
only on the binding energy of the dimer. 

In the dimer limit, a±/a^> 1 and 0# 3> 1, see Eq. ©, 
power counting shows that the first term in Eq. (|13fl is 
negligible. Using ({1/2, fl > 1) w -2%/0 and G(0,0) = 
4/3 yields from h(0) = — 1 the result 



dad = -KooV^sa± +/?aj_/i 
= — K 00 a]_/2a + 2/3a, 



(14) 



where Kqo = 9/32 = 0.28125. We also specify the sub- 
leading order in Eq. I|14|l. where (3 ~ 0.543. A similar 

calculation gives b a d/aj_ = (8/9)f2 B for 51^ 3> 1, vali- 
dating a repulsive zero-range ID atom-dimer potential in 
the low-energy limit, V ad (z) = g a d,S{z) with g ad > 0. Due 
to the presence of higher channels, a a d is renormalized in 
the dimer limit. This correction can be derived analyti- 
cally [lTj by making explicit contact to the integral equa- 
tion in the unconfined case [l3.ll5j. We find Hoc ~ 0.636, 
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FIG. 1: Scattering length a a d/a± versus dimensionless bind- 
ing energy fi_g. The solid curve is the numerical solution 
to Eq. ill .'-!> . and the dotted (dashed) curves represent the 
analytical results in the dimer (BCS) limit, respectively, see 
Eqs. ifHEjl and lfT5jl . 



which gives a a d = — (a^) 2 /2 (1.2 a) with the confinement 
scale a^ = (3h/2m uj±) 1 ^ 2 for the atom-dimer reduced 
mass 2mo/3. With the 3D atom-dimer scattering length 
1.2 a [13, this result exactly matches the dimer limit of 
the analogous two-body result J^J. Since closed channels 
do not cause profound changes even in the dimer limit, 
they are neglected in what follows. 

Outside the dimer limit, in general a numerical solution 
of Eq. H13(l is necessary. To ensure regularity of h(u), 
it is beneficial to Fourier transform back to real space, 
where the Fourier transformed h(z) is well behaved and 
allows for a quickly converging solution of the integral 
equation. The numerical result for a a d/a± as function 
of Ob is shown in Figure ^ I n the BCS limit, where 
a±/a <^i —1 and 0# <C 1, we find 

a a d = koci-l/v^b = Koa\/\a\, n w 0.75. (15) 

The solution in Fig. ^ for arbitrary O^ nicely matches 
onto the limits i|14|) and (|15fl . Remarkably, around 
Ob w 2.2, there is a zero of a a d- At first sight, this 
behavior seems to be linked to the two-body CIR, see 
Eq. J2J. However, the atom-dimer 'resonance' occurs at 
a different O^, and, more importantly, the assumption of 
a 5-potential interaction breaks down. This can be seen 
by computing the range parameter b a d m Eq. I|llf) . see 
Figure[3 While in the dimer limit, b a d stays small, in ac- 
cordance with our analytical result, in the BCS limit, it 

3/2 

is found to diverge as b a d oc O s . This implies that one 

cannot use a ^-potential to describe atom-dimer scatter- 
ing outside the dimer limit, but instead more complicated 
potentials have to be used, e.g., a repulsive square- well 
potential. Furthermore, the potential becomes non-local 
|19j . An effective ID potential for low-energy atom-dimer 
scattering can be constructed directly from Eq. ©. In 
fact, transforming Eq. JSJ back to real space, one gets a 
Schrodingcr equation 



— - P ) /(*) = - / dz'V ad {z, z')f(z') (16) 
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FIG. 2: Range parameter b a d/a±_ versus ils- The solid curve 
gives b a d from the numerical solution of Eq. 1131 . and the 
dashed line represents a a d from Fig. Q for comparison. 



with the non-local effective potential 



V ad {z, z) = - 



dk dk jkz-ik'z' k 
2ir 2ir 



C(Q k 



-A 



k,k' 



(17) 



where A^k' is given by Eq. JJjJ. It can be checked eas- 
ily that this potential becomes very wide in the BCS 
limit, with support given by the dimer size, but it always 
stays repulsive. Hence there is no three-body bound state 
(trimer) even in the confined problem for arbitrary a±/a. 
Finally, we have also analyzed the problem of three- 
body recombination and dissocation, T + T + I^TI + Tj 
see Ref. |l3j for the unconfined case, where a non-zero in- 
coming state ^o nas to be taken into account in Eq. 10. 
Antisymmetry of $o( x j) under y — > — y requires that 
the lowest transverse state (for y) has angular momen- 
tum m = ±1, implying that dissociation processes en- 
counter an energy barrier Hu>j_(l + fig) not present in 
the unconfined problem. The three-body recombination 



problem is therefore simpler, and dimers are rather stable 
against three-body dissociation. 

To conclude, we have solved the three-body problem 
of a binary cold Fermi gas. The scattering length aw 
describing atom-dimer scattering has been extracted for 
arbitrary confinement and 3D scattering length, and is a 
universal function of a±/a. For Qb ~ 2.2, correspond- 
ing to aj_/a w 2.6, we find aw = 0. However, this 
does not imply resonantly enhanced atom-dimer inter- 
actions. While in the dimer limit. a±/a ^$> 1, a standard 
repulsive zero-range potential is found, in the BCS limit, 
aj_/a <C — 1, the situation turns out to be more compli- 
cated. The scattering length is positive, aw > 0, but at 
the same time the range of the effective interaction be- 
comes very large, b a d 3> aw. Therefore it is not possible 
to employ zero-range potentials in that limit anymore. 
Nevertheless, it is worth pointing out that also the con- 
fined three-body problem is universal in the sense that it 
can be completely expressed in terms of two-body quan- 
tities. This is encouraging news for many-body calcula- 
tions which rely on model parameters extracted only from 
two-body physics [T^, D3 ■ Finally, we mention that the 
scattering length a a d (and possibly the length b a d) can be 
extracted experimentally using standard techniques, and 
thus the scenario put forward above could be checked in 
detail. Future work should also address the role of Pauli 
blocking due to the background Fermi sea. For the un- 
confined case, this was studied by Combescot J2(j. Fur- 
thermore, in a similar fashion as done here, dimer-dimer 
scattering in the confined geometry can be analyzed from 
the four-body problem. While in the deep dimer limit, 
contact to the free-space result [l(j can again be estab- 
lished, the behavior for Qb < 1 is particularly notewor- 
thy. This will be discussed elsewhere [l7j . 

This work was supported by the SFB TR12 of the 
DFG. 
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